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Decoherence in current induced forces: Application to adiabatic quantum motors. 
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Current induced forces are not only related with the discrete nature of electrons but also with its 
quantum character. It is natural then to wonder about the effect of decoherence. Here, we develop 
the theory of current induced forces including dephasing processes and we apply it to study adiabatic 
quantum motors (AQMs). The theory is based on Biittiker’s fictitious probe model which here is 
reformulated for this particular case. We prove that it accomplishes fluctuation-dissipation theorem. 

We also show that, in spite of decoherence, the total work performed by the current induced forces 
remains equal to the pumped charge per cycle times the voltage. We find that decoherence affects 
not only the current induced forces of the system but also its intrinsic friction and noise, modifying 
in a non trivial way the efficiency of AQMs. We apply the theory to study an AQM inspired by 
a classical peristaltic pump where we surprisingly find that decoherence can play a crucial role by 
triggering its operation. Our results can help to understand how environmentally induced dephasing 
affects the quantum behavior of nano-mechanical devices. 

PACS numbers: 85.85.-|-j, 73.23.-b, 03.65.Yz. 


I. INTRODUCTION. 


Nano-mechanical devices in general and nanomotors in 
particular are topics that have attracted much attention 
in recent years A— The working principle of most of these 
experimental and theoretical proposals rely on classical 
physics. However, at nanoscale one may benefit from 
the many phenomena emerging from quantum interfer¬ 
ences. In this direction, reverse quantum pumping has 
been recently proposed as the basic mechanism by which 
DC-current induced forces can drive a nanoscopic motor. 
Such a device is now known as “Adiabatic Quantum Mo¬ 
tor” (AQM)j2r— . Despite the classical description of the 
motor’s movement, one can profit from the quantum na¬ 
ture of current induced forces. Indeed, one can boost the 
efficiency of AQMs by exploiting the interferences present 
in a Thouless pumpZ. 

Many natural questions arise when one addresses cur¬ 
rent induced forces in nano-mechanical devices. Will 
their quantum behavior survive under realistic non-ideal 
situations? Which is the effect of environmentally in¬ 
duced decoherence and how can we model it? Will the 
effect of decoherence be always counter-productive for an 
AQM? Solving these issues will provide a better under¬ 
standing of their working principles and ways to assess 
the feasibility of their experimental implementations. 

In this work, the theory of current induced forces based 
on scattering matrices^Ai^ is extended to include deco- 
herent events. Our approach is based on a reformulation 
of Biittiker’s fictitious probe modeliA2&. This allows us 
to address the effect of decoherence on non-equilibrium 
current induced forces, friction coefficients, and fluctuat¬ 
ing forces with a focus on AQMs. 


II. THEORY. 

A. Langevin equation and cnrrent indnced forces. 

In the non-equilibrium Born-Oppenheimer approxima¬ 
tion, the dynamical degrees of freedom of a system are 
slow as compared to the electron dynamics. In this limit, 
we can treat the mechanical degrees of freedom as a clas¬ 
sical field acting on the electrons. Since, the movement 
of a motor is cyclic, one can reduce its many degrees of 
freedom to a single rotational coordinate, x. Hence, one 
can describe the rotor’s dynamics by the 1-D Langevin 
equation, 

Mx + ^=F-^x + ^, ( 1 ) 

dx 

where M is the mass of the rotor (or the moment of iner¬ 
tia), and U is some classical external potential that can 
also be introduced. The right hand side of Eq. [T] ac¬ 
counts for the current induced forces where F is a mean 
adiabatic reaction force. The second term is a friction 
(dissipative) force where 7 is the friction coefficient. The 
last term, accounts for the force fluctuations. These 
current induced forces have a quantum origin and prac¬ 
tical expressions in terms of the scattering matrices are 
given in RefsAiA^ which are our starting point. 

Like in Brownian motion, the interaction of the elec¬ 
trons with the rotor gives rise to a dissipation mechanism 
and a fluctuating force. At equilibrium, the only non¬ 
vanishing contribution to the friction coefficient is the 
symmetric contribution Qjj other hand, 

there is a fluctuating force f (t) whose self-correlation D, 
defined by (f(t)C(^O) ~ D6{t — t'), is assumed locally cor¬ 
related in time. These quantities must be related by the 
fluctuation-dissipation theorem (FDT), i.e., they satisfy 
D = 2KbT 7 ®’®?, where KbT is the thermal energy. We 
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will consider the forces only up first order in eV and/or x. 
That implies that only the equilibrium contributions to 
the friction coefficient and D will be taken into account. 


B. Decoherence and null current condition. 

Decoherence is included by connecting the system to a 
fictitious voltage probe (/). The voltmeter condition im¬ 
poses Irj) = 0, which provides charge conservation, while 
the reservoir character of the fictitious probe ensures 
the loss of memory of the re-injected electrons. Strictly 
speaking, this also involves inelastic events that redis¬ 
tribute electron’s energy. However, in a linear response, 
eV —>• 0, inelasticity is just reduced to a stochastic de¬ 
phasing of the wave function, i.e. decoherence.— 

The total current flowing at lead m is 

+ + ( 2 ) 


current cancellation of Eq. [2] is 

' u,+r„ y - —"v ■ 

( 5 ) 

Since the variation of x is slow respect to the electronic 
dynamics, we can consider that 5^^ adapts instanta¬ 
neously to satisfy /^ = 0 at all time. 

C. Current induced forces in presence of 
decoherence. 

By considering an infinitesimal bias, 6fj,L = eV, we can 
split the Fermi function into an equilibrium and an out- 
of-equilibrium contributions, fm = fo + ^fm- Then, the 
force generated by the electric current 



where is the non-equilibrium current caused by an 

infinitesimal bias eV, is the pumped current due 

to a variation of x, and 5Im accounts for current fluctua¬ 
tions. Fictitious volage probe metod has proved useful to 
address a quite related problem, the quantum pumping 
with dephasing— While there, fluctuations could be 
disregarded, here, we need to consider them explicitly as 
they affect the dynamics of the system. 

The non-equilibrium current incoming to the sys¬ 
tem through the lead m = L, in a linear re¬ 
sponse limit and at low temperatures, is 1^“* = 
(e/27r?i) 

m ■ Here, 

is the transmittance between leads m and n = L, R, (p] 
and Sfim = — is the chemical potential of m, taking 

fiR as reference. 

The pumped current through a lead m is 

jPump _ e(dnm/da;)i, were e is the electron charge, x 
the velocity of the rotational parameter x, and the emis- 
sivity of the lead m is 


dum 

dx 


de 

27ri 





( 3 ) 


where / is the Fermi distribution, H^ is a projector onto 
the lead m, and S is the scattering matrix— By inte¬ 
gration of we obtain the pumped charge per cycle 

through lead 


QL = e 


/ 


driL 

dx 


Tl^, dn^ 

Tl(I) + dx 


dx. 


( 4 ) 


Current fluctuations contain both the non-equilibrium 
shot-noise and the thermal noise. Both of them satisfy 
(dim) = 0, and ~ However, at equilibrium, 

the shot noise vanishes and only survives thermal noise. 

The null current condition /^ = 0 at lead (p directly 
imposes a condition to d^j,^. The value of d^j,^ that ensures 


can be split as F = + AF. The equilibrium force, 

is conservative and thus, it does not produce work. 
At low temperatures and for small eV, f {pAfmde ~ 
{■)dfj,m. Using Eqs. |3l [3 and M we calculate the non¬ 
conservative forces for systems without magnetic fields 
up to first order in eV and x. 


AF = 


driL Tl^ dn^ 

dx Tl 4 , + T^r dx 


eV 


—2Trh 


X 

Tl4, + T^r 



2Trh dl^ 

e -I- 


dn^ 

dx 


( 7 ) 


The first term on the right hand side (RHS) of Eq. [7] 
is the non-equilibrium force, F"®, whose second term 
within the brackets is the decoherent contribution. This 
force is responsible for the work of the system, which is 
obtained by evaluating W = j> {F^^ F F"®) dx, yielding 


W = 


dx 


driL 

dx 


TL(j> dn^ 
Tl 4 > + T^r dx 


eV. 


( 8 ) 


Comparing Eqs. [5] and |3] we realize that the work in 
presence of decoherence is proportional to the pumped 
charge per cycle times the voltage’s bias, 


W = QV. 


( 9 ) 


This relation, proved valid for coherent AQMs, can 
be thought as a signature of the Onsager’s reciprocal 
relations;^ showing that the model is well behaved. 

The second term on the RHS of Eq. [7] is a force pro¬ 
portional to the velocity. Thus, one can associate it with 
a dissipative force whose origin is purely due to decoher¬ 
ence. The resulting friction coefficient is 

^ _ 2'Kh f dn^ \^ 

Tl</> + t^r V dx y 


( 10 ) 
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FIG. 1. Scheme of the type of system studied. A rotational 
device interacting with a quantum dot. Rotation of the motor 
changes the dot’s energy as well as its coupling to one lead. 


This prove that decoherence inside the sample enables 
energy dissipation through an additional friction of the 

rotori22i2^ 

The third term on the RHS of Eq. [3 accounts for the 
fluctuations on the force induced by decoherence. The 
self-correlation of the fluctuating force can be defined 
as ~ D^6{t - t'). Thus, 




/ 2Trh 1 dn^ \ 
V e Tl^ + T^r dx ) 


( 11 ) 


where the spectrum power of the fluctuating current 
at the lead (f), is defined as ~ — In 

deducing Eq. 1111 we use that transmittances and emissiv- 
ities are already mean values. The only contribution to 
the current fluctuation which is non-vanishing at equi¬ 
librium is the thermal fluctuation or Nyquist-Johnson 
noised Thus, S,f, is characterized by 


= 2KbT ^ {Tl^ + T^r) . (12) 

ZlTfl 

Replacing Eq. [T^linto Eq. ITT] yields 

= 2KbT 7 -^, (13) 

which demonstrates that decoherent friction and fluctu¬ 
ating forces induced by decoherence are related by the 
FDT in our model as they should be. 


D. Efficiency of AQMs. 

The thermodynamic efficiency, Tym j can be defined as 
the useful output power that can be extracted from the 
system over the total input power. The useful output 
power is the work per cycle of the motor minus the en¬ 
ergy lost due to friction, all divided by the period r, i.e. 
QV/t—JJ "fx^dt/r. The total input power is IV+QV/t. 
Then, 


^TD 


Q - 47r^7*/ (rV) 
grV -f Q 


(14) 



on-resonance 


off-resonance 



FIG. 2. (Golor online) Upper Fig.- Total work of the AQM 
as function of (Eo — e) and F,^, both expressed in units of 
Fq. The work is normalized to its maximum value in the 
figure. Lower Figs.- Work of the AQM for the “on-” and 
“off-“ resonance conditions, (Eq — e) = 0 and —1 respectively. 
Different contributions are marked in the inset. In all figures, 
the work is normalized to the maximum value of the total 
work. 


Here, we have introduced I = gV where g is the average 
conductance, and the corrected average friction coeffi¬ 
cient 7 * = fi ^7 where 7 is the average friction coefficient 
7 = / -yx^dt/ j> x^dt. The dynamical constant d essen¬ 
tially accounts for the deviations of x respect to its mean 
value, d = t/{2tt) ^{x"^) where (x^) = j> x^dt/r. 

From Eq. [T4]one can extract the minimum energy nec¬ 
essary for rj ^ 0, which is the minimum energy necessary 
to start the motor’s motion QV > 47 r^ 7 */T. Addition¬ 
ally one can also realize that there is an optimal value of 
T that maximizes the efficiency. 


To 


47r^7* 

QV 



Q' 




(15) 


This value can be used to find the optimal load that the 
motor can move or, given a load, the optimal voltage to 
be applied to maximize the efficiency. Note that, when 
the average conductance goes to zero, tq goes to infin¬ 
ity, which is consistent with the adiabatic limit of the 
efficiency proposed in Ref?^. 
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FIG. 3. (Color online) Upper Fig.- Total friction coef¬ 
ficient normalized to its maximum value as function of Eq 
and both expressed in units of Fq. Lower Figs.- Fric¬ 
tion coefficient of the AQM for the “on-” and “off- “ resonance 
conditions, (Fq — e) = 0 and —1 respectively. Different con¬ 
tributions are marked in the inset. In all figures, the values of 
the friction coefficients are normalized to the maximum value 
of the total friction coefficient. 


III. RESULTS. 


Just to illustrate our theory, we will consider a simple 
example. Fig. [l] of a quantum dot that while rotating 
changes its resonant energy, E{9) — Eq +/S.E cos{9+ 9 q), 
and the coupling to one reservoir, Vr + AU sin(0). The 
angular coordinate is 9 and 9q is a phase shift. The cou¬ 
plings to the other reservoir L and to the fictitious probe 
(j) are assumed constant. The system-environment inter¬ 
action rate is determined by the coupling to the 

fictitious probe. The details of the solution of this ex¬ 
ample can be found in Appendix A. The main assump¬ 
tions used are: (1) The Landauer-Biittiker picture of non¬ 
interacting electrons is valid. (2) The interaction with the 
leads can be taken within the wide band limit (WBL). 
(3) The interaction between the dot and the movable 
part of the system is perturbative, i.e. the variations of 
AF and AU are small respect to Tq, the width of the 
dot’s resonance without decoherence at 0 = 0. (4) The 
terminal velocity of the rotor in the stationary regime is 
approximately constant, i.e. d ~ 1. We briefly discuss 
the conditions for the validity of this last in Appendix B. 

The work of the motor can be split into a coherent 
and a decoherent contributions. The behavior of W is 
similar to that described in Ref.^® for the (decoherent) 




FIG. 4. (Color online) Thermodynamic efficiency of the 
AQM shown in Fig. [T]as function of the decoherent rate F,^ 
(in units of Fq) and the period of the motor r (in units of 
fi/2ro). Upper Fig.- on-resonance regime. Lower Fig.- 
off-resonance regime (Fq — e = —O.TSFq). The efficiency is 
normalized to its maximun value in all plots. Insets show the 
cuts marked in the main figures as continuous or dot lines. 


pumped charge. Details of the solution can be found in 
Appendix A, but the main results are shown in Fig. [S] In 
the on-resonance regime, i.e. when the Fermi energy £ is 
e ~ Fq , the coherent contribution to the work is a mono¬ 
tonic decreasing function of F^, whereas the decoherent 
term increases with F ^ until it reaches a maximum value 
to decay afterward. The total work is always a decreasing 
function of F,^. In the off-resonance regime, the coherent 
and the decoherent contributions behave qualitatively as 
in the on-resonance regime. However, here the decoher¬ 
ent term dominates and, then, the total work presents a 
maximum for a finite value of F,^. This implies that in 
the off-resonance regime an adequate environment inter¬ 
action can indeed maximize the work of the motor. 

The friction coefficients and 7 “^ present a sim¬ 

ilar behavior both in the on-resonance and in the off- 
resonance regimes as function of F^, see Fig. [3l 7 ®’®'^ 
is a monotonically decreasing function of F,^ while 7 '^ 
presents a maximum. The total friction coefficient 
7 = 78.69 -(- 7 «^ always decays with F^, both in the on- 
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resonance and in the off-resonance regimes. In principle, 
just by looking Eqs. [7] and [13 one could (naively) ex¬ 
pect that electronic friction is increased due to the ex¬ 
tra friction term. However, electronic friction is actually 
mitigated by decoherence. This is because decoherence 
diminishes the quantum fluctuations and consequently 
also the friction coefficient, due to the FDT. 

The effect of decoherence on the thermodynamic effi¬ 
ciency is complex due to a competition between the total 
work per cycle and the friction coefficient. In fig. Ujwe 
plot the efficiency as function of T,^ and the period r. As 
predicted by Eq. |T3 there is a (decoherence dependent) 
optimal value of r that maximizes the efficiency. In this 
particular example the maximum efficiency is quite low, 
1.6 10“^. This is expectable considering the perturbative 
interaction between the rotor and the dot, which implies 
a small Q, and the high coupling between the dot and 
the reservoirs, which gives a high Beyond that, we 

find that, in the on-resonance regime, decoherence always 
diminishes the efficiency of the AQM. On the contrary, 
in the off-resonance regime, there is a maximum in the 
efficiency for a finite value of T^. This shows that, sur¬ 
prisingly, efficiency can be increased due to a properly 
tunned interaction with the environment. Indeed, there 
are certain values of r where at a critical T^ the effi¬ 
ciency switches from zero to finite values. This implies 
that the AQM is indeed triggered by its interaction with 
environment. 


IV. CONCLUSIONS. 


In conclusion, we developed the general theory of de- 
coherent current induced forces and we applied it to 
AQMs. We showed that decoherence not only modi¬ 
fies the current induced forces but also the electronic 
or intrinsic dissipation mechanisms and its related force 
fluctuations.— >2^ We proved that the theory is consistent 
with the fluctuation-dissipation theorem. Besides, we 
showed that the relation between the total work and the 
pumped charge per cycle, proved for coherent AQMs,— 
remains valid in presence of decoherence. We exempli¬ 
fied our theory with a simple example of AQMs showing 
that even there the role of decoherence can be non triv¬ 
ial. Indeed, we showed that decoherence can increase the 
efficiency of the AQM. This may allow the design of novel 
devices such as environmentally activated nanomotors. 
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Appendix A: Details of the calculation of the 
example of an AQM with decoherence. 


In this appendix we give the details of the calculations 
of the example presented in the main text, an AQM based 
on a peristaltic pump. The energy level of this AQM 
depends on 9, the rotational coordinate, as E{9) = Eq + 
AE cos{6 + 6 q), where 6q is a phase shift. The coupling 
to the reservoir R is Vr + AyRsin(0). We assumed a 
non-interacting electron picture. All the solutions are 
obtained considering a perturbative interaction between 
electrons and the mechanical degree of freedom, i.e. AE 
and AVr are small respect to Tq, the width of the dot’s 
resonance without decoherence at 9 = 0. The chemical 
potential of the right reservoir is taken as the reference 
energy and we assumed a wide band approximation. All 
observables are evaluated at low temperatures and in a 
linear response regime, where eV —>■ 0. All this implies 
£ = /ifl in all the equations. 

The effective Hamiltonian of the whole system, i.e., dot 
plus leads, reads 

H{9)=E{9) + El + ^r{9) + Y.^, (A1) 


where Ei is the self energy of the lead i = L,R, (j). We 
consider the leads within the wide band limit (WBL) 
where the self-energies E’s are pure imaginary quantities 
and are independent of the Fermi energy e. Thus, Er = 
-iPi, = -iP^, and Er{ 9) = -iT r{ 9) = -i[VR + 
AV/j sin(0)]^/Vg^, where 41/^^ is the bandwidth of the 
lead. Here, the quantities P r and P,^ are constants since 
the couplings to their corresponding leads are assumed 
independent of 9. The retarded Green function, defined 
as Go = [e — H{9)]~^, is 


e - iE{9) - iPi - iTR{9) - iT,^)' 


(A2) 


Transmittances Tm^n, needed to calculate the observables 
from Eqs. |l|to|3 can be obtained from 


|G„ 




(A3) 


where m ^ n, and m,n = L, R, </>. In this example, 
Gm,n = Go, for all m,n. We can also obtain the density 
of states (DOS), N{£,9) = —(l/7r)Im(Go), as 


N{£,9) 


1 m 

7r(e-A(6»))2-Hr(d)2’ 


(A4) 


where r(0) = Pl -I- T^ -|- ri^(d). The complete S-matrix 
can also be obtained from lA2l 

/ tr ^TRTRi9) yrzTV \ 

s = i- 2 iGo I yrlr^ Tr{ 9) y'TR{9)r^ . 

V ypITV ^TR{9)r^ / 

(A 5 ) 

Thus, the matrix [5’'!'^] gives 

5 t^ = - 2 i|Go|"A, (A 6 ) 
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where the elements of the operator A are 


Ai.i 

Ai ,2 


Ai,3 

A2,2 

^2,3 


A 3 ,3 


dE{0) 
^ d9 ’ 


dE{9) 

d9 


{e - E{9) + ir(0)) 


^/r^ 


4rK(0) 
dE{9) 


d9 


d9 




d9 


d9 




P AE{9) 

^ d0 ’ 


dTR{9) 

d9 


with Ajn,n = 

At this point, we can evaluate all the physical quanti¬ 
ties that are relevant to this problem. Let us start with 
the emissivities of Eq. [31 At low temperatures, we can 
use dfmide = -/m(l - Imj/KsT ~ -d{e - ^im)- Thus, 
we have 


dn^ 

d9 




1 

27ri 



where m = L, R, 4>. Then 


dni(0) 

d9 

duR 

d9 


N{e,9))dE{9) 
r(6») d9 
l N{e,9)) 

TT r(6») 

dE(6») 
d6» 


^Rid) 


L{4>)^ 


drfl(^) 

d6i 


ie-Ei9)) 


(A7) 

(A 8 ) 


(A9) 


We can insert these expressions into Eq. 0] to obtain the 
pumped charge per cycle or the total work, note that 
W = QV. The result can be split into a coherent and 
a decoherent contributions, W = W‘^°^ + The 

mathematical expression for is exactly the same 

with or without decoherence (however, its value depends 
on r^). Consistently, is zero in a purely coherent 

case. 

By evaluating Eq. [8] with Eq. IA9I and using Green’s 
theorem;^ we obtain 


W = -eVi^coh + fldec)A.EAVR cos(6»o), 


where 




^dec — 


TrV;^ 
T^TlVr 

ToETEg^ 


27rAf(£:) 


To 


+ 4Tr^N'^{e) 


(AlO) 

(All) 

(A12) 


Here, Tp = r(0)|AVii=o,A£;=o which does not depend on 
9, and Tg = ^These equations are expressed 
at zero order, i.e., for small AE and AVr. Note that the 
factor AEAVr cos(6*o) is the parametric area enclosed in 
a cycle of the parameter 9. 

The friction coefficient at equilibrium, jg 

m,n 

xTr|n^5t^n„5t^|. (A13) 

Assuming a low temperature limit to evaluate the inte¬ 
gral, yields 



The additional friction term, which is a direct conse¬ 
quence of decoherence on current induced forces, can be 
evaluated from Eq. (TUI giving 

7^ = 

Notice that positivity of both coefficients are guaranteed. 

List of parameters. The parameters used in the whole 
work to perform the calculations are: eV = 10“"^ 

Er = VOT Eg^, Tr = 0.1 Eg^, AEr = 10-2 Vr, 
AE = 2 10-3 Eg,^, and Tg = Tr + Ej/Eg^ = 0.2 Eg„. 
The maximum values of the efficiency, total work, and 
total friction coefficient used to normalize the plots are 
1.6 10-3, lO-'^eV, and 2.4 10-3?i respectively. 


Appendix B: Limit of constant terminal velocity. 


At steady state, all the energy that is absorbed by the 
motor is completely dissipated by friction and thus, total 
energy is conserved according to 


Fi9) - Eioad - 


d9 = 0. 


(Bl) 


where Fioad account for an external load to the system. 
Here, we are assuming an average over steady state en¬ 
sembles, so random forces can be neglected and we can 
use 9{t) — 9{t -I- r) = 6{9). 

From eq. El] and using the same definitions for 7 and 
(d‘^{t)\ as in the main text, we obtain. 


^02 (t)^ = (B2) 

where Wioad is the work done by the external forces in 
one cycle. Using this result, the mean kinetic energy 
becomes 

{K) = 1/ (9\t)) = . (b3) 
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where I is moment of inertia of the rotor. 


netic energy of the motor is much greater than At/* {9 ), 


l^ (QV-W^load) 

2 7 


» At/*(0), 


(B4) 


The change of the kinetic energy due to a rota¬ 
tion of the parameter 9 is AK{9) = —AU*{9), where 

U*{9) = - \f{9') - Fioad - j0{9')] d9' is a pseudo¬ 

potential defined only at the steady state. When the ki- 


then AK{9)/ (K) —?> 0 and, thus, the rotational veloc¬ 
ity of the system becomes insensitive to 9. Therefore, 
T = d9/9{9) ft! 27r/9 and hence the dynamical factor 

~ Note that the terminal velocity is 

independent of /. Therefore, the condition given in Eq. 
IB4I are always valid in the limit of macroscopic rotors, 
rotors with a sufficiently large moment of inertia. 
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